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Abstract. Let T„ be the nilpotent group of real nxn upper-triangular matri- 
ces with Is on the diagonal. The hamiltonian flow of a left-invariant hamilton- 
ian on T*T„ naturally reduces to the Euler flow on t* , the dual of in = Lie(r„). 
This paper shows that the Euler flows of the standard riemannian and sub- 
riemannian structures of T4 have transverse homoclinic points on all regular 
coadjoint orbits. As a corollary, left-invariant riemannian metrics with posi- 
tive topological entropy are constructed on all quotients D\T„ where D is a 
discrete subgroup of T„ and n > 4. 



1. Introduction 

Let E be a nilmanifold, i.e. homogeneous space of a connected nilpotent Lie group 
G. Each homogeneous riemannian metric on G induces a locahy-homogeneous 
metric on E. These riemannian geometries, which wih be called left-invariant, are 
of interest in both geometry and dynamics. A basis question is 

Question A: Which left-invariant geodesic flows on a compact nilmanifold have 
zero topological entropy? 

A mistaken answer to question A appears in Theorem 3 of [9j. In j2], the first 
author showed that on 2-step nilmanifolds, all left-invariant geodesic flows have 
zero entropy. In 0| , metrics on compact quotients of the 3-step nilpotent Lie group 
Ti © T3 are constructed whose geodesic flows have positive topological entropy. 
The paper also speculated that the standard geodesic flow on T4 also had such 
horseshoes. Montgomery, Shapiro and Stolin [10] investigated the standard suh- 
riemannian geodesic flow on T4 ^ ; they showed that it reduces to the Yang-Mills 
hamiltonian flow which is known to be algebraically non-integrable [141 115] . 

Let us state the first result of the present paper. The Lie algebra of T4, {4, has 
the standard basis consisting of those 4x4 matrices Xij with a unit in the i-th 
row and j-th column, i < j, and zeros everywhere else. A quadratic hamiltonian 
/i : 14 — > R is diagonal if it is expressed as h{p) — J2i<j '^ijiPj -^ij)'^ fo'" some 
constants aij. The standard riemannian metric has aij — 1 for all i, j; the standard 
Carnot (subriemannian) metric has ai2 = 023 = 034 ~ 1 and all other coefficients 
zero. 

Theorem 1.1. If h : R is a diagonal hamiltonian with 012013023034 7^ and 

O13O34 ~ 012024, then for all but at most countably many regular coadjoint orbits in 
i\, the Euler vector field of h has a horseshoe. In particular, the Euler vector field 
of the standard riemannian metric (resp. sub -riemannian metric with 013 ^ 0) is 
analytically non-integrable. 

The condition that 013O34 = 012O24 is only a device to simplify the proof: all 
nearby hamiltonians also have a horseshoe. We also show, by means of a numerical 
computation of an integral (see Table [T] in section Hj), that when 013 = 0, the 
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conclusions of Theorem 11.11 hold. This shows that the subriemannian geodesic 
flow of Montgomery, Shapiro and Stolin is real-analytically non-integrable. Related 
numerical computations (figure [1] in section |4]) suggest that the hamiltonian h has 
a horseshoe on every regular coadjoint orbit provided only that 012023034 ^ and 
013034 = 012O24. Let us formulate a corollary to Theorem 11.11 Let -D < T„ be a 
discrete subgroup of T„, S = -D\T„ and ST. is the unit sphere bundle. 

Theorem 1.2. If n> A, then there is a left-invariant geodesic flow <f>t ■ ST ST 
such that htop{4>i) > 0. 

It appears likely that all left-invariant geodesic flows on ST have positive topo- 
logical entropy and are non-integrable with smooth integrals. 

Theorem 11.21 is interesting from a riemannian point of view. Let {M, g) be a 
smooth (C°°) riemannian manifold, and ipt ■ SM SM the geodesic flow of g. 
For each T > and p,q G M, let nT{p, q) denote the number of distinct geodesies of 
length no more than T that join p to q. Mane [S] showed that if AI is compact then 
htop{(f'i\SM) = limT^oo ^os/mxa/ ^t{p, q) dpdq. Thus, for the geodesic flows 
constructed here, for generic points p and g on a compact quotient of T„, nxip, q) 
grows exponentially fast. In constrast, Karidi showed that the volume growth on 
the universal cover of these manifolds is polynomial of degree ^n{n^ — 1) [S]. 

Theorem 11.11 is proved by reducing the hamiltonian flow of (/>t on T*T to a 
hamiltonian flow on the coadjoint orbits of Lie(r4)*. In the appropriate coordinate 
system, the reduced hamiltonian is a small perturbation of a hamiltonian on R** 
that is the sum of an unforced Duffing hamiltonian and a forced linear system whose 
solutions can be expressed in terms of Legendre functions. The Poincare-Melnikov 
technique developed in [71 lllj for autonomous Hamiltonian systems is adapted here 
to show that a suspended Smale horseshoe appears in the perturbed hamiltonian 
flow. 

2. The construction on LiE(r4)* 

In this section, we will first recall a number of key facts about geodesic flows 
and left-invariant hamiltonian systems on the cotangent bundle of a Lie group; for 
more details, see [S]. We will then reduce the equations of motion of a left-invariant 
geodesic flow on T*T4 to the equations of motion of a hamiltonian system on T*R^. 

2.1. Poisson geometry of left-invariant hamiltonians. A Poisson manifold 
is a smooth manifold M such that C°°{M) is equipped with a skew-symmetric 
bracket {, } that makes (C°°(M), {, }) into a Lie algebra of derivations of C°°(M). 
The centre of (C°°(Af), {, }) is traditionally called the set of Casimirs. If / is a 
Casimir then Xf = Q and / is a first integral of all hamiltonian vector fields. If 
the set of Casimirs of (C°°(Af), {, }) are the constant functions, then we say that 
(C°°(M), {, }) is a symplectic manifold. In this case, the Poisson bracket naturally 
induces a closed, non-degenerate skew 2-form on M which is called a symplectic 
structure. We will say that a smooth map / : M ^ TV is a Poisson map if 
/* : {C°°{N),{,}n) -> (C°°(M), {,}m) is a Lie algebra homomorphism. 

The most basic example of a Poisson manifold that is also symplectic is provided 
by r*R — {{a,A) : a,A G R} equipped with the Poisson bracket satisfying 
{a,A}T'R = 1. 

The dual space of a Lie algebra gives an example of a Poisson manifold that 
is not (in general) a symplectic manifold. Let g be a finite-dimensional real Lie 
algebra and let g* be the dual vector space of g. T*g* is identified with g for all 
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peg*. The Poisson bracket on g* is defined for all f,he C°°(g*) and p € fl* by 
(1) {f,h}{p):=-{p,[dfp,dhp]), 

where (•, •) : g* x g ^ R is the natural pairing. Recall that for ^ G g, ad^ : g* ^ g* 
is the linear map defined by (ad^p, ry) = — (p, [^,77]). ad* : g gl{Q*) is the 
representation contragredient to the adjoint representation. For any h € C°°{q*), 
the hamiltonian vector field Eh = {•, /i} equals — ad^^^p. The standard example of 
a hamiltonian vector field is obtained from a positive-definite linear map (p '■ Q* ^ Q 
by setting h{p) = l{p,(j){p)), in which case Eh{p) = -ad^(p)p. 

Let G be a connected Lie group whose Lie algebra is g. The adjoint represen- 
tation of G on g, Adg^ = ^|t=o5exp(i^)5~^, induces the coadjoint representation 
{Ad*p,^) = {p,Adg-i£_) for all p & q* , g & G and ^ G g. As each vector field 
p — !■ ad^p is hamiltonian on g*. with linear hamiltonian h^{p) = —{p,£,), the coad- 
joint action of G on g* preserves the Poisson bracket. The orbits of the coadjoint 
action are called the coadjoint orbits. Each coadjoint orbit is a homogeneous G- 
spacc, and every hamiltonian vector field on g* is tangent to each coadjoint orbit. 
For this reason, the Poisson bracket {, }g. restricts to each coadjoint orbit, and is 
non-degenerate on each coadjoint orbit. Thus, the coadjoint orbits are naturally 
symplectic manifolds. A Casimir is necessarily constant on each coadjoint orbit, 
and in many cases (as in this paper) each coadjoint orbit is the common level set 
of all Casimirs. 

The Poisson bracket on g* also arises in a natural way from the Poisson bracket on 
T*G. The group G acts from the left on T*G, and this action preserves the Poisson 
structure {, }t*g- The set of smooth left-invariant functions C°° {T* G)^ is therefore 
a Lie subalgebra of C°°{T*G) with respect to {, }t*g- This subalgebra is naturally 
identified with (C°°(g*), {, as follows: the left-trivialization of T*G = G x g* 
induces the projection map r : T*G — > g* onto the second factor; r*G°°(g*) = 
G°°{T*G)'^ and r* is a Lie algebra monomorphism. 

The hamiltonian flow of a left-invariant hamiltonian H on T*G therefore has the 
equations of motion: 

g = TeLgdh{p), 

p = -a.d*ihip) P, 



(2) XH{g,p) 



Note that dr(Xij) = Eh, where h € G°°(g*) satisfies r*h = H. The vector field 
Eh is called the Euler vector field. It is a natural reduction of Xh by G. If 
h{p) = for a positive-definite linear map : g* — * g then H is induced 

by a left-invariant metric on T*G and Xh is the geodesic vector field. 

Finally, if £> < G is a discrete subgroup, then T*E = E x g* where S = D\G. 
The projection map r : T*G ^ g* is naturally left-invariant, so it factors through 
to a map Vo : T*E g*. If Ho = r*h for some h e G~(g*) then TroiXnJ = Eh- 
Thus, the hamiltonian flow of a left-invariant hamiltonian on T*E always projects 
to a hamiltonian flow on g*. 

2.2. Poisson geometry of T*T4. Let denote the Lie algebra of T4, so 

X z w 

y u 



u 



V 




u, V, w,x,y,z gIL 



For each a e {u, v, w, x, y, z}, let A e t4 be the element obtained by setting a equal 
to one and all other coefficients equal to zero. Then {U,V,W, X,Y, Z} is a basis 
of t4 whose commutation relations given by: [AT, F] = Z,[Y,V] = U, [X,U] = W, 
[Z, V] = W, and all others are trivial or obtained by skew-symmetry. 
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Let Pa : fl* ^ R be the linear function given by Pa{p) — {p, A) for A E g and all 
p G 2*. From the definition of the Poisson bracket on t| (c./. Eq. [l]), along with the 
commutation relations, we conclude that: 

{Px,Py} = -Pz, {Pv,Pv} = -Pu, {Px,Pu} = -Pw, {Pz,Pv} = ~Pw 

There are two independent Casimirs of 14 are Ki{p) — p^, K2{p) = PwPy—PzPu- 
Let K : tl ^ be defined by K = {Ki, K2). The level sets of K are the coadjoint 
orbits of T4S action on and will be denoted by Ok, where k = (fci, fe). We will 
say that Ok is a regular coadjoint orbit if fcifc2 ^ 0. 

Lemma 2.1. Let Ok be a regular coadjoint orbit. Then Ok is symplectomorphic 
to (T*R)^ equipped with its canonical symplectic structure. 

Proof. The Poisson bracket on t| restricts to Ok. We will denote the restricted 
bracket by {,}fe. Let (r*R)2 = {{a,A,b,B) : a,A,b,B e R}. The canonical 
Poisson bracket, which we denote by [, ], satisfies [a, A] = 1, [b, B] — 1 and all other 
brackets are zero. 

Let A and fi be two non-zero parameters (the parameters are included because 
we will further transform coordinate systems). The map fk{p) = (a, A, 6, B) defined 

by 

a = -Xpx, A = {kiXy^pu, b = -fipv, B = (fci/i)"Vz 
is a diffeomorphism of Ok onto r*R^. Indeed, fk is clearly smooth. And gk{a, A, b, B) — 
{pu, ■■■,Pz) defined by 

Pv = -IJ-^^b, pu, = ki, px = -A^^a, py = (/c2 + klXfiAB)/ki, and p^ = ki^B, 

satisfies K o gk = k and fk o gk = id, gk o fk = id. Since gk is an algebraic map, it 
is smooth, and so we see Ok is diffconiorphic to T*R^. 

The commutation relations for the Poisson bracket on allow us to compute 
that {/fefl, fkA}k — {fkb, fk^}k = 1 and all other Poisson brackets are zero. It 
follows that fl : (C°°(r*R2), [,]) ^ (C°°(Ofc), {, }k) is a Lie algebra isomorphism. 
Hence fk : Ok T*I{? is a symplectomorphism. □ 

2.3. The hamiltonians. Let > be constants such that ai3034 = ai2a24 and 
let 

(3) 4_ff (p) = ai2p^ + a23P^ + aiap^ + a24P?. + a34p^ + auP^ 

Since the vector field Eh is unaffected by the addition of a Casimir, the term ai4p^ 
can be ignored. 

Let us introduce a symplectic change of variables on T*R^: A = -^{X — Y), 
B — -^{X + Y) , a = '^i^~y)' b = "71 ^ ^'^^ Z = C. Because 013034 = 

012024, there exists unique X, ^ > so that = 034/1^^ — oi2A^^ — 013/1^ — 024A^, 
and oi2A^^ + 034/x^^ = 1. Indeed, we can choose A^ = 2oi2 and /i^ — 2034. Then: 

(4) 2Hfe = {x^ - ^X^ + vX^) + + uoY"^ + vY^ - 2vX'^Y^), 

where ^ -(013034^1+023^:2^012034), lo = ^+2013034/0^ = ai3034fcj-023fc2 ^'^12034 
and V — 012023034/;;^. Note that we can write u = ^ + 2cv where c = . 

Lemma 2.2. The Euler vector field of H on the regular coadjoint orbit Ok (equation 
is a time change of the hamiltonian vector field of 

(5) 2H = + (^X^ - + a^Y^ +Y^ - 2X^Y^ 
on r*R^, where = 1 + 2cj/3 . 
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Proof. Define the coordinate change g,y{x, X,y,Y) = [ax, a ^X,ay,a ^Y) where 
a = i/s. Then g*(Hfe) = a^H. □ 

Lemma 2.3. For all e > 0, the hamiltonian flow o/H (equation\Bj) is conjugate to 
the flow of the vector field 

X = X, Y ^ y, 

X = X -2X^ + 2eXY'^, y = [-a^ + 2X'^] Y + 2eY^ . 

Proof. Introduce the coordinate transformation h^{x, X,y,Y) = (x, X, y^y, y^Y). 

□ 

Remark. It is clear from ^ that the vector field depends on the parameter a. 
Inspection of the formula for a (lemma l2.2l and immediately above) shows that a is 
identically unity when the coefhcient ai3 vanishes. This is the case for the standard 
subriemannian metric, where ai2 — 023 — a^^ — 1 and the other coefficients vanish. 
Rather than specializing to a = 1, we have elected to carry a through our analysis. 
The rationale for this will be apparent in section [3.31 



(6) Xe = 



3. Analysis of 

For e = 0, the map is singular. However, the vector field Xq is well-defined. 
We will show that Xq has a normally hyperbolic invariant manifold S whose stable 
and unstable manifolds coincide, and that this manifold S persists for e > but 
the stable and unstable manifolds W^{S) no longer coincide. This implies that the 
Euler vector field EnlOk has transverse homoclinic points for all regular coadjoint 
orbits. 

3.1. The normally hyperbolic manifold S. Inspection of the vector field X^ 
shows that the set 

S = {{x,X,y,Y) : x = X = 0}, 
is invariant for all e. One sees that for e = 0, the vector field is 



X = X, Y ^ y, 

X = X + 0{X^), y = [-a^ + 2X^]Y, 



which shows that S is normally hyperbolic. Therefore S is normally hyperbolic for 
all £ sufficiently small. Since X^ is conjugate to the same vector field for all non-zero 
e, one concludes that S* is a normally hyperbolic manifold for all e. 



2 _ 1 \2 
2) 



3.2. The stable and unstable manifolds of S. The function h = x'^ + {X 
is a first integral of Xq. The set h~^{j) is the stable and unstable manifold of S, 
which we denote by W^{S). On W^{S) - S, the flow of Xo satisfies 

X = ± sech(i -I- to), x — =F tanh(t + <o) sech(t -I- io)^: 

Y = coYQit + to) + ciYi{t + to), y ^ Y, 



(7) 



where X{0) = ± sech(io), x{0) — ^ tanh(io) sech(io)^ and Yj solves the initial- value 
problem 



(8) 



f + - 2sech(i)2] r = 0, (*) 
r(0) = l-j, Y{0) = J 



while y(0) = CQYo{to) + ciYiito), y{0) = coio(to) + ciYi{to). The solutions Yj are 
chosen so that they are even (j ~ 0) and odd (j = 1) functions of time. 
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3.3. The Melnikov function. To determine if the flow of has transverse ho- 
moclinic points for non-zero e, we appeal to the following theorem. 

Theorem 3.1. Let ip^ : M x R — > ill be a complete, smooth flow that depends 
smoothly on e. Assume that (po possesses a normally hyperbolic, invariant manifold 
S C M , and that there is a smooth function h : M R such that 

(1) the stable and unstable manifolds of S coincide and equal h^^{0); 

(2) dh does not vanish on W^{S) — S. 

Then, for all sufficiently small non-zero e, p^ possesses a normally hyperbolic in- 
variant manifold Se and the local stable and unstable manifolds of (W^{S) and 
W~{S), respectively) can be written as the graph of a function sf : W^{S) — > 
W^{S). The splitting distance, defined for p g W^{S) by s^{p) — ho sf{p) ^ ho 
s~{p), is a smooth function of e and s^(j>) — em(j>) + O(e^) where 

(9) m(p)= / {dh,y)opl{p)dt, 

and y = — — 
oe at 



Proof. The proof of this theorem is a standard application of invariant manifold 
theory plus an adaptation of the proof of the Melnikov formula [3 [TTJ |3] . □ 

Remark. If m changes sign, then, for all e 7^ sufficiently small, the perturbed 
stable and unstable manifolds intersect but do not coincide. In our case, a result 
of Burns and Weiss implies that the topological entropy is non-zero. If m has a 
non-degenerate zero, then the implicit function theorem implies that, for all e 
sufficiently small, W^{S) has a transverse intersection with W~{S). Note that 
the intersection of the surface S with a constant energy level is a periodic orbit. 
Therefore each trajectory in the intersection is doubly asymptotic to a periodic 
orbit in S. 



For the flow defined by X^, we have that 

X = 0, y = 0, 

X = 2Xr2, y = 2Y^. 



(10) y 



Whence 

(11) {dh,y) = AxXY"^ , 

since h — x^ + [X — . The equations in ([7|) imply that the Melnikov function is 

(12) m{p) = cqCi X / -4 tanh(T) sech(T)^ Yo(7-) Yi(t) dr. 

Remarks. (1) The formula for the Melnikov integral (fT^ appears to be a function 
on S not Wq(S) — S. This does not contradict Theorem (|3.ip . Inspection of the 
integral (equation [^1) shows that m((^g(p)) = m{p) for all s and p. The coordinate 
system on W^{S) — S determined by ([7]) uses time along the flow as one coordinate 
{to), so only the other two coordinates, cq and ci, ought to appear in the Melnikov 
function. (2) If we write m{p) = 2coCi x /, where 

/•CO 

(13) /= / -4 tanh(T) sech(T)2 ro(T) Yi(t) dr, 

then m has non-degenerate zeros along {cq = 0, ci 7^ or ci — 0, cq 7^ 0}, provided 
that / 7^ 0. 
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3.4. The Legendre functions and /. Substitution of z = tanh(t) transforms the 
differential equation ^) into the Legendre differential equation 

(14) (1 - z^)Y" - 2zY' + (^^.(i. + 1) - -J^^ Y = 0, 

where /i — ia, v = + -^^^ and ' — j^- The integral / (equation [T3|) is trans- 
formed to 

(15) / = / zUoiz)Uiiz)dz, 

Jo 

where Uj{z) = Yj{t). 

3.5. The Melnikov function is non-zero: 1^0. For the remainder of this note, 
q : [0, oo) — s- R is a continuously differentiable function such that limt_,oo qit) = 
and a > is a fixed positive number. The function z = z(t) is assumed to solve 

(16) z+[a^ - q{t)] z = 0. 



In analogy with the integral / in equation p3p . define an integral 

/•OC 

(17) /= / g(t)zo(t)zi(t)dt 



/ is implicitly a function of a; one wants to prove that / can vanish at most 
countably many times. 

Let us first prove the following. 



Lemma 3.2. If z solves equation U6]) . then z is bounded with bounded derivative. 

Proof. Define H — ^{a^z'^ + P). From (fTH]) one computes that H — qzz. Inte- 
grating by parts yields H = Co + ^q{t)z{t)^ + /J —^q{s) z{s)^ ds, where Co is a 
constant that depends only on z(0) and z{0). Thus 

(18) i [a^ _ q(t)] z{tf <Co + J^ -\ m z{sf As. 

Since q{t) ^ as i — > oo, there is a T > such that — q{t) > ^a^ for all t > T. 
Therefore, equation implies that there is a constant Ci such that for alH > 

(19) z(t)2 < Ci + 4j X / -q{s) z(s)2 ds 

Jo 

This is a Gronwall inequality for u — z'^. Thus 

z{t)' < Ci exp{~^q{t)) 

for all t > 0. Since q is continuous and converges to at infinity, it is bounded. 
Therefore, z is bounded. 

To prove that w — z is bounded, define H = ^{a^uP' + uP'). Since w + 
\o? — q[t)\w — f, f — qz, one computes that H ~ qww + fw. One can bound 
Jq fis)w{s)ds using that z — w and that q and z are already bounded. One then 
obtains a Gronwall inequality like p9p for w(t)'^. □ 

Lemma 3.3. If Zo,zi are solutions to equation il6]) . then the limit 

(20) / = - lim \zo{t) zi{t) + a^zoit) zi{t)] 

exists and equals Wacot{B) where the angle B is defined below, and W is the 
Wronskian of the solutions zq, zi. 
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Proof. Let z be any solution of p^ . Let i„ be the sequence of zeros of z{t), indexed 
in increasing order. Let 0„ — atn mod 2tt and a„ = z{tn). There is a sequence 
such that 0„j, ~^ (f) mod 27r and a„j, ^ a > 0. The former foUows by compactness 
of R/27rZ and the latter because i is bounded. 

Since q{t) — > 0, the Sturm comparison theorem implies that the sequence 
converges to and t„+i — t„ converges to 7r/a [13] . One sees that positive and 
negative zeros of z must alternate for all n. Hence, without loss of generality, one 
may assume that z{t2n) — > a and i(i2n+i) — > —o.- The continuous dependence of 
solutions on initial data therefore implies that z{t + converges in the 

weak Whitney topology to a sin ( at— 0). In particular, z"(t) converges uniformly 
to a sm{at — 0) for t G [0, Air /a] . 

To apply these observations to the limit ([20|) . let N be sufficiently large so that 
the ri-th and n + 1-th zeros of both zq and zi are at most ^ apart for all n > N. 

Let s e [— , ^'"+^^^ 1 and write s = t + so that t e [0, Then 

|-io(s)ii(s) + 0:^2:0 (5)2:1(5) — a^aofli cos((?!)i — (j)o)\ 
< jig (t)i"(t) — a^aofli cos(af — ^o) cos(ai — 0i)| + 
a'^|2Q (t)2"(i) — flofli sin(Q;t — (po) sin(Q;t — (f>i)\- 

If s — » c», then n 00. The above-mentioned uniform convergence for t G [0, 27r/a] 
shows that the limit ([20]) exists and equals ^cos(-B) where A = a^aoai and B = 
01 — 00 mod 27r. 

On the other hand, the Wronskian W oi zq, z\ is constant and 

W = ZQ{t)zx{t) - zi{t)zo{t) ^ aoofli sin(0i - 0o), 

by the same argument as above. Therefore I/W — acot{B). □ 



Remarks. (1) The angle B has the following interpretation which emerges from 
the proof of lemma p.3p . The zeros of solutions to ([TO]) are asymptotically ir/a 
apart, and the zeros of linearly independent solutions are interlaced. The angle 
B is defined so that B/a mod 7r/a is asymptotically the time between consecutive 
zeros of the linearly independent solutions. Figure [U left, plots i? as a function 
of a for the solutions zj = Yj to the initial- value problem ([S]). One expects that 
as a ^ 00, the solutions should converge quite quickly to cos and sin, whence 
B should approach ^. The figure captures this behaviour quite nicely. (2) The 
function / — W^acot(i?) from lemma [3731 can be computed numerically. The Sturm 
comparison theorem implies that the n-th zero i„ of a solution z satisfies n/a < 
tn+i — tn < Tr/a X (1 + q/a'^) if \q{t)\ < for all t > t„. If q goes to zero 

sufficiently fast, one can numerically compute the first several zeros and obtain 
a reasonably accurate estimate of B. Figure [H right, shows the graph of / for 
q = 2sech(t)2. 

Lemma 3.4. // zq,zi are solutions to 116\) such that zq and z\ vanish at t = 0, 
then the integral 

/•oo 

([73) / = / q{t)zo{t)zi{t)dt 

Jo 

exists and equals Wacot(B) where the angle B is described in Lemma \3.S\. and W 
is the Wronskian of the solutions 29,21. 
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Proof. By lemma [XH both solutions are bounded, so one can apply integration by 
parts to the integral. This yields 

POO 

I = g(0)zo(0)zi(0) - / q{t) [zo{t) zi{t) + zo{t) h{t) ] dt, 

JQ 

q{t) [io(<)zi(<) +Zo(t)ii(t)] dt 

since zi vanishes at t = 0. 

From (|16p . it is known that q(i)zo(i) = za{t) + a^zo(i) and similarly for zi. 
Therefore 

f °° d 

I = - J — [zo{t)zi{t) +a^ZQ{t)zi{t)] dt, 

— — hm \zo{t) zi{t) + a'^zo{t) zi{t)] , since io(0) = = 2:1(0), 
= Vl^acot(B) by lemma 

□ 



Lemma 3.5. Assume that there exists C, A > such that \q{t)\ < Ce'^* for all 
t > 0. Then the integral I = /(a) is a holomorphic function of a on the strip 
|Ima| < A about the real line. 

Consequently, if q is an even, monotone function, then I vanishes countably 
many times at most. 

Proof. A solution z = z(t; a) to p6p is a holomorphic function of a for each fixed 
t [T^. For large t and |Ima| < A, the solution z ~ z(t;a) is equal to cos(ai + (p) 
plus a term that grows slower than e"^*. This implies, by the residue formula, that 
I = I {a) is holomorphic provided that jlmal < A. 

When a — and q is even, the even and odd solutions to (fT6| do not change 
sign. Therefore, if q is monotone, then the integrand defining /(O) does not change 
sign, so /(O) ^ 0. Thus, / can vanish at most countably many times on the strip 
|Ima|<A. □ 

Theorem \l.l\ If ai3 7^ - whence c 7^ in Lemma 12.21 -. then lemma shows 
that the hamiltonian flow of H (equation [3]) on all but countably many coadjoint 
orbits in i\ has a horseshoe. This proves the main result, Theorem ll.il □ 

4. The degenerate case when a = 1 

If ai3 = 0, as occurs for the Carnot subriemannian metric of |10j . then a = 1 
and lemma cannot be applied. We investigate two distinct ways to address this 
problem. The first is direct and numerical; the second leads to some further insight 
into the integral /. 

4.1. Numerical evidence. In this case, figure [T] indicates that /(I) is approxi- 
mately —2.75. Table [1] shows the results of a numerical computation of /(I) with 
varying step sizes. It is clear from this table that /(I) = —2.76 to two decimal 
places. 

To estimate the error in the computations, one uses the fact that the differential 
equations ()16p are hamiltonian with the hamiltonian 

(21) H = ip2_^i[a2-g(r)] z^ + u, 

where p, z and u, t are canonically conjugate variables (along solutions, r = tq + t, 
so it is a pseudo-time). Since H is preserved along solutions to (jl6p . the maximum 
deviation of H along a numerical solution provides an estimate of the upper bound 
of the error in the solutions zq,zi. 
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4.2. Qualitative evidence. As explained in the Remark in subsection 13.51 one 
may compute / as a function of a by computing the pliase angle B. Figure [T] 
graphs B and / versus a. This figure shows that /(I) does not vanish. 




Figure 1. Left: Phase angle B vs. a; Right: I — acot{B) vs. a. 
Both plots use the solutions Zj = Yj oi ^ with q{t) — 2sech(t)^. 
These solutions are computed numerically in Maple using the 4-th 
order Runge-Kutta method; the zeros are located by interval halv- 
ing. 
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Table 1. The numerical calculation of / with a = 1. 



0.5 

0.25 

0.125 

0.0625 

0.03125 

0.015625 

0.0078125 



'O.max 



Hi. 



mm 



Hi. 



max 



Ho.r: 



Hl,r^ 



Hi.,: 



-2.76812630 
-2.76366763 
-2.76340793 
-2.76339200 
-2.76339101 
-2.76339095 
-2.76339094 



-0.5 
-0.5 
-0.5 
-0.5 
-0.5 
-0.5 
-0.5 



-0.49025150 
-0.49944022 
-0.49996571 
-0.49999786 
-0.49999986 
-0.49999999 
-0.49999999 



0.49833857 
0.49992738 
0.49999582 
0.49999974 
0.49999998 
0.49999999 
0.49999999 



0.51067514 
0.50063022 
0.50003878 
0.50000241 
0.50000015 
0.50000000 
0.50000000 



0.00974849 
0.00055977 

3.42847126 Xl0~^ 

2.13207876 xlO"'' 

1.33088170 xlO-^ 

8.31541421 xlO"" 

5.19672801 xlQ-l" 



0.01233657 
0.00070283 

4.29572646 xlO^^ 

2.67223559 xlO"" 

1.66835298 xlO"'' 

1.04238692 xlO~** 

6.51456639 XIO^IO 



O 

CO 

< 

i 

O 

O 
d 



Solutions to the hamiltonian equations of H arc computed with the Forest-Ruth 4-th order symplectic 
integrator [T^ and initial conditions z(0) — j, i(0) = 1 — J, t(0) = 0, w(0) = for j ~ 0, 1. The maximum 
(resp. minimum) value of H along the j-th numerical solution over the interval [0, 35] is indicated by 



H, 



(resp. Hj 



The integral / is computed by 



N 



4=0 



where z'' is the computed solution with step size h, N = 35//i, ti = i x h and q{t) = 2sech(t)^. Files at 



http : / / www ■ maths . ed . ac. uk/~lbutler/t4 . html 



CO 

O 

•2, 
I — [ 

g 
> 

o 
d 



